Local Operations in qubit arrays via global but periodic Manipulation 
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We provide a scheme for quantum computation in lattice systems via global but periodic manip- 
ulation, in which only effective periodic magnetic fields and global nearest neighbor interaction are 
required. All operations in our scheme are attainable in optical lattice or solid state systems. We 
also investigate universal quantum operations and quantum simulation in 2 dimensional lattice. We 
find global manipulations are superior in simulating some nontrivial many body Hamiltonians. 
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I. INTRODUCTION 

Since a quantum computer (QC) will exhibit advan- 
tages over its classical counterpart only when a large 
number of qubits can be manipulated coherently, many 
architectures of QCs based on scalable physical sys- 
tems have been investigated 0, 0, S S 13- Among 
these candidates, ultra-cold neutral atoms trapped in the 
periodic potentials of an optical lattice are attracting 
much attention!! H H H U, 03 OH E2 Since neu- 
tral atoms couple weakly to the environment, decoher- 
ence is suppressed greatly. Besides, some typical many- 
body models can be constructed in optical lattice sys- 
tems. This makes it possible to use quantum optical 
methods to study fundamental condensed-matter physics 
problems [H[ll El [H. 

In the schemes for QCs based on optical lattice sys- 
tem, the couplings between nearest neighbor atoms in 
the lattice can be simultaneously switched on and off 
by adjusting the intensity, frequency, and polarization of 
the trapping light. However, it is difficult to focus a laser 
beam on a single atom due to the short lattice period. 
This makes it challenging to realize controlled collisions 
and perform single and two qubit logical operations. A 
possible approach to overcome the difficulty of addressing 
the atoms individually is based on marker qubits 0,0,0, 
i.e., by using two types of atoms with different internal 
states in the optical lattice. One type of atoms act as 
marker qubits to address logical qubits. Another method 
is one way quantum computation[l7| . In this proposal, 
the couplings between adjacent atoms can be utilized to 
prepare the atoms in a high dimensional cluster state. Af- 
ter the initialization, one can enlarge the lattice period to 
single qubit addressable range. Universal quantum com- 
putation can then be implemented simply via a series of 
single qubit measurements. Of course, one way quantum 
computation requires a great qubit overhead. 

Recently, a series of schemes for quantum computation 
based on global operations only were proposed. In these 
schemes, local addressability become unnecessary [13. Il9l 
I20L I2H I22I |23| . Remarkably, Raussendorf devised a 
novel scheme for universal quantum computation via 
translation-invariant operations on a chain of aubits[2lj|. 
in which only the nearest neighbor Ising-type interaction 
and translation-invariant single qubit unitary operations 



are required. Inspired by these proposals, in this letter 
we present a scheme for quantum computation based on 
global operations in qubit arrays only. Here, we develop a 
general method to realize single qubit unitary operations 
via global but periodic manipulation. Compared with 
the quantum cellular automata scheme [l8l Il9| . encod- 
ing overhead is eliminated in our proposal. In addition, 
single qubit operation is easier than in Raussendorf's 
scheme |21|. 

The paper is organized as follows. In Sec. II, we show 
that universal quantum computation in one dimensional 
arrays can be realized via global but periodic manipu- 
lation. We also display that the preparation for initial 
states and the measurement for final states can be imple- 
mented by using this global manipulation. Furthermore, 
in Sec. Ill, we discuss quantum computation and quan- 
tum simulation in high dimensional arrays. We find the 
implementation of the periodic operations in the finite 
lattice will benefit from global manipulation. Section IV 
contains some conclusions. 



II. UNIVERSAL QUANTUM COMPUTATION 
IN ONE DIMENSIONAL ARRAYS 

Let us first consider a two-component bosonic atomic 
mixture trapped in one dimension by two spin-dependent 
lattice potentials. Each atom is assumed to have two rel- 
evant internal states, which are denoted with the effective 
spin index a =f, J,, respectively. In the Mott insulator 
regime, the system can be described by a two-component 
Bose-Hubbard Hamiltonian||: 

H = - y^ftp-qt q i+lff + H.c.) + i ^ U a n ia {n ia - 1) 
+U}i^n l] n ii . (1) 

i 

Here, af a and are creation and annihilation opera- 
tors for bosonic atoms of spin a localized on-site i, and 
n-ia = a\ a ai a . Under the condition that t <C U a ,U\i and 
( n it) + ( n ii) = h Eq. (1) is equivalent to the following 
XX Z type interaction Hamiltonian|5l l2^| : 

H * = E [ J i«+i - + aVa* +l )] . (2) 
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Here, erf = n iT - n a , of = <4j.<z;x + afja lT , and 
- -i(a^a 4 | — a^flit) satisfy erf, erf 
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j! = (if + tf)/2c% - tf/c/ T - q/u u J 2 = t^tju u . 

Remarkably, the parameters J\ and J 2 can be easily con- 
trolled by adjusting the intensity of the trapping laser 
beam or an external fieldQ- Therefore, the following 
well-known spin 1/2 interaction models can be realized: 
Ising model (Ji ^ 0, J 2 = 0), XY model (Ji = 0, J 2 ^ 0), 
and Heisenberg antiferromagnctic (ferromagnetic) model 
(Ji = ±J 2 ). 

Besides Hamiltonian Hs, in order to implement uni- 
versal quantum computation, we introduce the following 
effective periodic magnetic field: 



2m 

COS I h tp 



Here, we set the distance between the nearest neighbor 
sites as unity, a is the period of the magnetic field, when 
a — ► 00, Hg(a,tp) reduces to a homogeneous magnetic 
field. To realize this periodic magnetic field, we set the 
eigenstates ||) and ||) to have the same energy. By the 
left and the right circularly polarized light, they are sep- 
arately coupled to the common exited level |e) with a 
blue detuning A. In such a 3-level system, we can ob- 
tain the effective Hamiltonians a x , a y , and er 2 in the 2 
dimensional Hilbert space spanned by |f) and ||) only 
by adjusting the polarization of coupling light HJ. Thus, 
for the cold atoms trapped in one dimensional lattice, we 
may construct such a periodic potential field Hg(a,tp) 
along the lattice direction by applying a monochromatic 
standing wave laser beam. In addition, in our scheme, 
we require the ability to adjust the period of the effective 
magnetic field Hg(a,ip). To this end, we split a beam 
of monochromatic light into two beams and make them 
interfere. Thus, we can obtain a one dimensional stand- 
ing wave along the interior bisector of the angle between 
two beams of lights, whose period depends on the angle 
between the two beams of lights. 



A. Single qubit operations 

Here, we present a general method to control any 
single qubit on one-dimensional finite lattice by using 
global Hamiltonian Hg(a,(p). Our idea for local ma- 
nipulation derives from Fourier transformation i.e., any 
spatial function can be decomposed into a superposi- 
tion of periodic eigen-functions. Let us consider a one- 
dimensional chain of N qubits. The unitary evolu- 
tion controlled by the effective periodic magnetic filed 
H%(a,<p) is: Ui{0,a,k,a) = exp {-if H %{a, 
Our aim is to implement a delta function type of oper- 
ation on the finite spin chain by assembling these wave- 
like unitary operations Ui(9,a,k,a). Without loss of 
generality, we consider performing a unitary operation 
cxp(— iOaf,) on the fcth site. The operation on a single 



site can be realized in an iterative way. In the iterative 
procedure, a key step is to construct parity unitary oper- 
ation: = U 1 (ir/A,2,k,z)Fj; 0) Ui(n/4,2,k,z)Fj; 0) = 
Y[\ k ~ 1 ]<j<\ N - k ^ exp(—i26af. 1 2j -). Here, [77] refers to 

the integer part of the number 77 and = 
exp {— i-^Hg(oo, 7r)}. The parity unitary operation 

F^ only remains homogeneous on even sites ar- 
ray from the fcth site. It effectively just imple- 
ments a homogeneous operation on a spin chain of 
N/2 qubits. We find that similar parity-eliminated 

operations can be established iteratively: F^' = 

U 1 (7r/A,2 l ,k,z)F^- 1) Ul(TT/4,2 l ,k 1 z)F i k l ~ 1 \ 1 < I < 
m,m = |~log 2 (N — k)~\ +1. Finally, operation on the fcth 

qubit can be realized: F^ = exp(— i2 m 9a% ) (see Fig. 
1(a)). In this iterative procedure, implementing single 
qubit rotation takes 3 x 2 m — 2 elementary steps (appli- 
cations of Ui(6, a, fc, a)), which roughly ranges from 1.57V 
to 67V. 
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Figl (a) 



Fig 1 (b) 




Fig 1 (c) 



FIG. 1: Schematic for addressing and readout measurements 
via global operations, (a) The process of single qubit op- 
eration. In an iterative step, approximately half sites are 
selected out. (b) A high efficiency procedure for operation 
between two nearest neighbor sites. In an iterative step, ap- 
proximately two thirds couplings in the remaining spin chain 
are eliminated, (c) Mapping quantum states |1) localized in 
far apart sites into auxiliary energy level \ai) by using modi- 
fied periodic magnetic fields. The read out for state \ai) can 
be implemented by the method of detecting the fluorescence. 



B. Two-qubit operations 

Since single qubit operations can be realized, im- 
plementing two-qubit operation between two near- 
est neighbor qubits via many body interaction be- 
comes feasible. We consider the case of Ising model: 

Hising = Z^i 1 J i a i a i+i- When the couplings be- 
tween nearest neighbor qubits switch on, the cor- 
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responding unitary evolution is Ui{6 = Jit) = 
Y\i exp(i#<T?<7?, j). To implement the interaction be- 
tween the kth and (k + l)th qubits we may synthesize 
the following unitary transformations: exp(i9a k a k+1 ) = 

^(!K^(-!K^i^(-jK^(!K4 + r 

In the above manipulation, a two-qubit operation takes 
6 single qubit operations, so 97V <~ 36 N elementary 
steps (applications of Ui(-d,a,k,a) or Ui("3)) are re- 
quired. This overhead can be greatly reduced by de- 
vising a new focusing procedure, which is shown in Fig. 
1(b). This focusing procedure is a bit similar to that 
of single qubit operation. The first iterative step is: 

T { k\e) = f/i(f , 3, k-1, x)Ui(0)U;(%,3,k-l, x)Ui{9) = 

■ ■ ■ exp(i29a z k _ 3 a^ 2 ) exp(i2ea z k a z k+1 ) exp(i2ea z k+ ^a z k+4 ) ■ 

■ ■ . Here, the script k refers to focusing the interac- 
tion between the kth and (k + l)th qubits. After the 
first step, almost two thirds couplings in this Ising chain 
are cut off. Then, we construct the iterative procedure: 



T { k n \9) = U 1 {a n ,3 n ,k-^,x)T^- 1 \e)U^(a n ,3 n ,k- 

^-,a;)T^ 1) (0). Here, a n is determined by the follow- 
ing restrictions: 



3 «-i + 1 x /3"" 1 -! 

COS j — 7T 1 — COS ( — — 

3" - 1 



2p+l 



COS 



3™ 



7T — COS 



3™ 



I 3"-l 



2 
m 



7T, (4) 



where p, I, and m are arbitrary integers. Obviously, for 
any given precision, we can always find an a n to sat- 
isfy Eq. (4). Finally, after repeating the iterative steps 
rif (<~ 0(log 3 N j) times, which corresponds to 3 x 2 nf — 2 
elementary steps, we obtain our anticipated operation: 
exp(i2"/0<7^ +1 ). 



C. Initialization and read out measurements 

To prepare the qubit trapped in each lattice site to |0), 
we drive the system to Mott insulator regime with one 
atom occupancy per lattice site. Furthermore, we adjust 
the interaction Hamiltonian Hs to Ising type. When 
Ji < 0, Ising type Hamiltonian have twofold degener- 
ate ground states |00,...,0) and |11,...,1). By applying 
a homogeneous magnetic field J2if a i (/ < 0)> the de- 
generacy of ground state will be broken. Eventually, the 
system will relax to ground state 1 00, 0) if the environ- 
ment is cold enough. 

To realize read out measurements for the final state 
of quantum computation, we may take advantage of 
auxiliary energy levels of the atoms (see Fig. 1(c)). 
Our strategy is to transfer internal state |1) to aux- 
iliary energy level |a;) only for qubits localized in far 
apart sites, which can be achieved by using modi- 
fied focusing operations. Therefore, we may address 
and read out internal states |a/) in these lattice sites. 



Compared with single qubit operation, we replace pe- 
riodic magnetic field J2i A [l — cos + <p)] erf with 
£\ A [l - cos (^et + <p)] <?ki, where a ki is Pauli opera- 
tor in Hilbert space spanned by {\ak~i) i , (we set 
l a o)i = Thus, we may map |1) to \a\) for atoms in 
next nearest neighbor lattice sites after the first iterative 
step. Using similar iterative procedure, we may transfer 
|1) to \ai) only for atoms 2 l sites apart. For the internal 
state | a;), we may address and read out by the method 
of detecting the fluorescence. 



III. QUANTUM COMPUTATION AND 
QUANTUM SIMULATION IN HIGH 
DIMENSIONAL ARRAYS 

Since focusing any single qubit operation and Ising in- 
teraction between any two nearest neighbor qubits can be 
realized, universal quantum computation can be imple- 
mented in one dimensional qubit array via global opera- 
tions. In the described scheme, a quantum computation 
on N qubits roughly takes O(N) (exactly less than 6 AT) 
elementary operations compared with addressable quan- 
tum computation schemes. Although an exponentially 
increasing number of steps can be avoided, a linear de- 
pendence on the number of qubits is still a great obstacle 
for large-scale quantum computation due to the fragile 
many body coherence. To reduce the required resources 
further we may consider quantum computation in high 
dimensional periodic qubit arrays. 

Actually, focusing processes in high dimensional 
lattice are quite easy because they can be realized just 
by using one dimensional focusing methods, repetitively 
Here, we just investigate quantum computation on 
two dimensional square lattice. In principle, higher 
dimensional quantum computation can be easily realized 
in a similar way. We set the number of sites in the row 
and column to n, Thus, the total number of sites in this 
lattice is n 2 . For simplicity, we denote the site of the 
ith row and jth column by First, let us consider 

a single site on two dimensional lattice. Without loss of 
generality, we give the unitary operation exp(— i#cr? .j) 

on site X(i,-%,x)Y(j,z,z)X(i,%,x)Y(j,-^z), 



where X(i,^,x) 



i"lE=i". 



nLi ex p(- i l c7 fi,fe)) and 

In this process, only 



0(n) elementary steps are needed. Next, we 
consider implementing Ising interaction between 
sites and + Since focusing operation 

on a single site is attainable, a straightforward 



way to realize the operation exp(z6*crf, , ) a i 



(i,j) u (i,j+i) 



^^(-iKjj^j+D^dKjj^j+i,- Here, 
U'j(9) = e X p{-iej2k<iW z {k ,i)^k±i,i) + ^*,i) CT (fc,i±i)]} 
In this focusing procedure, 8 operations on single 
sites are required. To reduce the operational re- 
sources, we can also consider the simplified scheme. 
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As shown by the above section, we can use the more 
efficient proposal to implement the interactions be- 
tween the nearest neighbor sites in one dimension, i.e. 
to produce an operation exp{— i0^2k a (k j) a (k i+i)^ 
via 0(log 3 n) iterative steps. Hence, we just 
need 2 single site operations to achieve our goal: 



From the above proposal, we may find that the quan- 
tity of operations will be saved further for quantum com- 
putation in high dimensional qubits arrays. Similar anal- 
ysis will show that, compared with addressable quantum 
computation schemes, 0( d ^/N) times numbers of oper- 
ations are required for the scheme in a d-dimensional 
qubits array via global operations. Here, N refers to the 
total number of the qubits. Then, we may make a rough 
comparison for local addressable and global addressable 
schemes from the viewpoint of the error of the unitary 
evolution. We set the single step average error rate e 
for local operations and £ for global ones. For simplicity, 
we look on e and £ as the upper bounds of two types 
of single step errors. As far as the scheme via global 
operations is concerned, an overhead of 0( d y/~N) times 
numbers of operations is cost compared with address- 
able ones. While, Ref[25[ shows that the upper bound 
of the error of the unitary evolution will linearly increase 
with the number of unitary transformations. Therefore, 
if e > 0( d VN)^ one can benifit from the scheme via 
global operations. While, in the above analysis, we do 
not differentiate the upper bounds for the different types 
of operations, which will lead to more subtle comparison 
for the overhead of operations in some well-known algo- 
rithms. A further analysis for the upper bound of the 
error is necessary and meaningful. In addition, an inter- 
esting question remains open: what is the best strategy 
for quantum computation via global operations? 

Once universal quantum computation can be imple- 
mented, any SU (2™) unitary transformation can be 
achieved. Therefore, in principle, quantum computers 
have the ability to simulate dynamical behavior of any 
finite dimensional quantum systems. However, it is a 
formidable task to decompose the time-dependent uni- 
tary evolution under some non-trivial many body Hamil- 
tonians into a sequence of elementary quantum gate op- 
erations. Fortunately, one can use Trotter formula to 
approximately implement quantum simulation: 



iEU^t = lim TT, 

n — ► no l A A 



-iHit/n 



(5) 



\i=l 



But, it is still quite difficult for a standard quantum 
computer to simulate some highly correlated many body 
models. As an example, let us consider "coupled dimer" 
Hamiltonian 1261 : 



H d = J 



XJ £ 



(6) 



■*■■*■■*■■* 
■*■■♦■■*■■* l=> 
■*■■*■■*■■* 
■*■■*■■*■■* 



Fig 2 



FIG. 2: The coupled dimer antiferromagnetic Hamiltonian 
Hd is represented by the two dimensional lattice on the right. 
The A links are shown as thick solid lines, and the B links 
as thick dashed lines. Hd can be decomposed into the square 
lattice antiferromagnetic part and decoupled dimer part. 



Here, A links form decoupled dimers while B links 
couple the dimers (see Fig. 2). A is the dimensionless 
coupling coefficient. To simulate the unitary evolution 
exp{— iHdt} of this 2-dimensional lattice with n 2 sites, 
it will take 0{n 2K ) elementary logical gate operations, 
if we partition [0, t] into K identical intervals. However, 
the overhead can be avoided if we use global operations. 
We may decompose Hamiltonian H d into two parts: 



H d = H 1 +H 2 , where H 1 = XJJ2 



H 9 



(i-A)JE 



and 



Hi is the square lattice 



(».i>e a 



b 



anti-ferromagnetic Hamiltonian, which can be directly 
realized by adjusting the potential of the optical lattice. 
In addition, it is quite straightforward to implement uni- 
tary transformations U a (9) = exp{i9 j\^A a i a j} an< ^ 
U b {6) = cxp{/rtV , ; ; ,;T;rr: + afa v 3 )} by using hori- 
zontal global operations and parity unitary operations. 
Because of ctJ commutes with afa^+a\a v -, unitary trans- 
formation exp{^ jjSjyj Hz} can be achieved by concate- 
nating U a (9) and [/& (9). Hence, more effective simula- 
tion for many body Hamiltonian H d can be implemented 
by using global manipulation. By investigating this ex- 
ample, we find it is very convenient to use global manipu- 
lations to implement operations with periodic structures, 
although focusing processes on single or two sites will 
take many iterative steps. This suggests that quantum 
error-correction can be achieved by taking advantage of 
global manipulations. 

In addition, by adjusting the period and phase of 
the Hamiltonian described by Eq. (3) we can imple- 
ment the unitary operation Ilfe=o hkhk+icrlk^^lk+s- 
Thus, an interval interaction form can be realized: 
Ylk=o ex P(^ (T 2fc cr 2/c+i)- A s shown by Ref[l9l|. if one can 
realize four types of Hamiltonian: -ff^li+i' 
and H^Li 2i in one dimensional lattice, we can imple- 
ment quantum computation via cellular-automata ap- 
proach. Here, #| i(2i+1) and fl^2i +1(M _ li2i) refer to the 
uniform single qubit operation on even (odd) sites and 
the interaction between the even sites and its right (left) 
nearest neighbor sites, separately. 
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IV. SUMMARY 

In this paper, we have devised a scheme for universal 
quantum computation in the periodic lattice via global 
but periodic operations. Our idea for quantum manipu- 
lation derive from mathematical Fourier transformations, 
i.e., any a function with space distribution can be decom- 
posed into the superposition of periodic eigen- functions. 
As a result, we derive the controls localizing single site 
and two nearest neighboring sites by adjusting global 
periodic Hamiltonians. We also show that simulating 
nontrivial many body Hamiltonian and quantum error- 



correcting operations will benefit from global manipula- 
tions. 

Note added: After submission of this work, we noted 
two schemes for implementing single-qubit operations in 
2 dimensional optical lattice [271 128| . 
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doctoral Fellowships, Z-W. Zhou thanks L.-M. Duan, L.- 
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